Lineshape predictions via Bethe ansatz for the one-dimensional spin-1/2 Heisenberg 

antiferromagnet in a magnetic field 
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The spin fluctuations parallel to the external magnetic field in the ground state of the one- 
dimensional (ID) s = i Heisenberg antiferromagnet are dominated by a two-parameter set of 
collective excitations. In a cyclic chain of N sites and magnetization < M z < N/2, the ground 
state, which contains 2M Z spinons, is reconfigured as the physical vacuum for a different species 
of quasi-particles, identifiable in the framework of the coordinate Bethe ansatz by characteristic 
configurations of Bethe quantum numbers. The dynamically dominant excitations are found to 
be scattering states of two such quasi-particles. For N — > oo, these collective excitations form 
a continuum in (g,tj)-space with an incommensurate soft mode. Their matrix elements in the 
dynamic spin structure factor S zz {q,u}) are calculated directly from the Bethe wave functions 
for finite N. The resulting lineshape predictions for N — ► oo complement the exact results 
previously derived via algebraic analysis for the exact 2-spinon part of S zz (q,u}) in the zero- 
field limit. They are directly relevant for the interpretation of neutron scattering data measured 
in nonzero field on quasi- ID antiferromagnetic compounds. 
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I. INTRODUCTION 

Advances in experimental techniques combined with 
improvements in sample preparation make it possible to 
produce data of ever increasing resolution for the quan- 
tum fluctuations and the underlying collective excitations 
in quasi-one-dimensional (ID) magnetic compounds. Ad- 
vances in the theoretical analysis of relevant model sys- 
tems combined with progress in the computational treat- 
ment of aspects that remain elusive to exact analysis 
make it possible to gain an ever more profound under- 
standing of the observable collective excitations in terms 
of a small number of constituent quasi-particles. 

There is scarcely a better case for illustrating this 
multi-track advancement of understanding quantum fluc- 
tuations than the ID s = | Heisenberg antiferromag- 
net and the growing number of materials that have been 
discovered to be physical realizations of this model sys- 
tem. The Hamiltonian for N spins h arranged in a cyclic 
chain with isotropic exchange coupling J between nearest 
neighbors and a uniform magnetic field h, 
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H — ^ [JS„ • S n+ i — hSn] 



(1) 



is amenable to exact analysis via Bethe ansatzQi and 
displays dynamical properties of intriguing complexity. 
The field h is a controllable continuous parameter, which 
leaves the eigenvectors unaltered, but changes the nature 
of the ground state via level crossings and thus has a 
strong impact on the dynamical properties, in particular 
at low temperatures. 



At h > hs = 2J the ground state of H has all spins 
aligned in field direction: \F) = | f j • • • j) is the reference 
state of the coordinate Bethe ansatz, and all eigenstates 
are described as excitations of interacting magnons, a 
species of spin-1 quasi-particles. Hence \F) plays the role 
of the magnon vacuum. The ground state \A) of H at 
h = contains N/2 magnons. The Bethe ansatz en- 
ables us to reconfigure this state as the physical vacuum 
for a different species of quasi-particles - the spinons, 
which have spin |. The entire spectrum of the Heisen- 
berg model (|l|) can_also be generated as composites of 
interacting spinons n 

Both descriptions are valid throughout the spectrum, 
but the magnon interpretation is more useful near the 
magnon vacuum, and the spinon picture is more use- 
ful near the spinon vacuum. The interaction energy of 
magnon scattering states or spinon scattering states is of 
0(N^ 1 ) as long as the number jof quasi-particles in the 
collective excitations is of 0(1 ).□ In a macroscopic sys- 
tem, the spectrum of such states is thus indistinguishable 
from the corresponding free quasi-particle states. Even 
under these simplifying circumstances, however, the in- 
teraction of the quasi-particles remains important in the 
make-up of collective wave functions, and is likely to 
strongly affect the transition rates and lineshapes. 

At intermediate values < h < hs of the magnetic 
field, the number of magnons or spinons contained in 
the ground state of H is of 0(N), implying that the in- 
teraction energy for either quasi-particle species remains 
nonzero for N — > oo in the ground state and in all low- 
lying excitations. This obscures the role of individual 
magnons or spinons in the collective excitations and ob- 



2 



structs the interpretation of spectral data obtained by 
experimental or computational probes. 

We can circumvent this problem by configuring the 
ground state \G) at < M z /N < \ as the physical vac- 
uum for yet a different species of quasi-particles. From 
the new vantage point, the dynamically relevant collec- 
tive excitations are then again scattering states of few 
quasi-particles with an interaction energy of 0(iV~ 1 ), 
which greatly facilitates the interpretation of the spec- 
tra probed experimentally or computationally. 



II. DYNAMIC STRUCTURE FACTOR 

In an inelastic neutron scattering experiment per- 
formed at low temperature, the observable scattering 
events predominantly involve transitions from the ground 
state to a subset of collective excitations filtered from the 
rest by selection rules and transition rates. Under ide- 
alized circumstances, the scattering cross section is pro- 
portional to the T = dynamic spin structure factor 



S^(q,cj) = 27rY / \(G\S^\X)\ 2 S(u;~ 

\ 



(2) 



where S£ = A^ 1/2 e l9n S*£, M — x iVi z IS the spin 
fluctuation operator. In a macroscopic system, the ag- 
gregate of spectral lines in (||) pertaining to scattering 
events with energy transfer u>\ = E\ — Eq, momentum 
transfer q = k\ — kc, and transition rate |(G|5^|A)| 2 form 
characteristic patterns of spectral weight in (q, w)-space. 
The shape of the spectral weight distribution provides 
key information on how the dynamically relevant col- 
lective excitations are composed of quasi-particles with 
specific energy-momentum relations. 

Experimentally it is possible, at least in principle, to 
separate the information contained in the dynamic struc- 
ture factors of the spin components parallel and perpen- 
dicular to the field direction, i.e. the functions S zz (q,u)) 
and S xx (q,u>) = \[S+-(q,w) + S- + (q,w)], respectively, 
for the fluctuation operators S% and S 1 ^ = ± iS]j. 
At h = the additional symmetry of H dictates that 
S zz {q,u)) = = ±S-+(q,w). 

An anchor point for the new results presented in the 
following is the exact 2-spinon dynamic spin structure 
factor at T = 0, which was determined recently via al- 
gebraic analysis and shown to contribute 73% of the to- 
tal intensity in S zz Lg, uj) at h = OE Given the energy- 
momentum rclatior 



- Jsinp, < p < 7r, 



(3) 



of the spinon quasi-particle, the 2-spinon states with 
wave numbers q = p\ + p 2 and energy oj — e S p(pi) 
e S p{P2) form a continuum confined by the boundaries! 



e Uq) = ^ J l sin, 2i eu(q)=TT.J 
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FIG. 1: Exact 2-spinon lineshapes at q = ir/2,3ir/4,ir 
of S zz (q,u>) at T = for Hamiltonian (Eh at h = as 
determined via algebraic analysis. The inset shows the 
boundaries (f|) of the 2-spinon spectrum. 

as illustrated in Fig. [l] (inset). The main plot shows the 
exact 2-spinon lineshapes of S zz (q, uS) at q = 7r/2, 37r/4, tt. 
The most detailed experimentaljdata available for testing 
these results pertain to KCUF3.B 

We shall see that the magnetic field causes dramatic 
changes in both the spectrum and the lineshapes. At the 
root of these changes is a change in the nature of the rele- 
vant quasi-particles. Two compounds suitable for study- 
ing magnetic-field effects on spectrum and lineshanea are 
Cu(C 6 D 5 COO) 2 -3D 2 and Cu(C 4 H 4 N 2 )(N03)2.&N 



III. BETHE ANSATZ EQUATIONS 

The Bethe ansatzB is an exact method for the calcu- 
lation of eigenvectors of integrable quantum many-body 
systems. The Bethe wave function of any eigenstate of 
(Q) in the invariant subspace with r — N/2 — M z reversed 
spins relative to the magnon vacuum, 



\ip) = a(n x , 

l<ni<... <n r <N 

has coefficients of the form 



,n r )S ni ■■■S nr \F), (5) 



a(ni, ... ,n r ) 



ves 



exp 




ViVj 



(6) 



(4) 



determined by r magnon momenta ki and one phase angle 
6ij = —Oji for each magnon pair. The sum V G S r is over 
the permutations of the labels {1,2,... , r}. 

The consistency requirements for the coefficients 
a(ni,... ,n r ) inferred from the eigenvalue equation 
H\ip) — E\ip) and the requirements imposed by transla- 
tional invariance lead to a set of coupled nonlinear equa- 
tions for the fej and A computationally convenient 
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rendition of the Bcthe ansatz equations has the form 
N<l>(zi) = 2irI i + '22<l>[(zi-z j )/2], i = l,...,r, (7) 

where 4>(z) = 2arctanz, ki = it — (fi(zi) and Oij — 
7rsgn[3fi(zi — Zj)] — (f)[(zi — Zj)/2\. Every solution of 
(0) is specified by a set of Bethe quantum numbers 
I\ < I2 < • • • < I T 5 which assume integer values for 
odd r and half-integer values for even r. The energy and 
wave number of the eigenvector thus determined are 



E-E f 
J 



St 



i=l 



(8) 



where Ep — JN /A is the energy of the magnon vacuum. 

We consider the class K r of eigenstates whose Bethe 
quantum numbers comprise, for < r < N/2 and < 
m < N/2 — r, all configurations 

~ + \-m<h<I 2 <--- <lr<\-\ + ra. (9) 

Here we employ the solutions {zi} of the Bethe ansatz 
equations not only to generate spectral data via (||), 
which is standard practice, but also to evaluate tran- 
sition rates |(G|5^|A)| 2 for the dynamic structure fac- 
tor (0) directly from the normalized Bethe wave func- 
tions |A) = The compiilational aspects of this 
method are discussed elsewhere. Lil 



IV. PHYSICAL VACUUM AND 
QUASI-PARTICLES 

The ground-state wave function \G) at < M z < N/2 
is specified by the set of r = N/2 — M z Bcthe quantum 
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(10) 



As the magnetic field increases from h — to hs = 2J, 
the magnetization M z increases in units of one from zero 
to N/2. A sequence of level crossings produces a magne- 
tization curve (M z /N versus h) in the form of a staircase 
with N/2 steps of height-Jj/iU- which converges toward a 
smooth line as N -> oo-BBEJ 

Depending on the reference state used for the charac- 
terization of the ground state \G), it can be regarded as 
a scattering state of N/2 — M z magnons excited from 
the magnon vacuum \F) or as a scattering state of 2M Z 
spinons excited from the spinon vacuum \A). To illus- 
trate the distinct roles played by the two species of quasi- 
particles in the class- K r states, we show in Fig. || the con- 
figuration of Bethe quantum numbers for \G) in a system 
with N = 8 and all values of M z realized between h = 
and h — fig. The positions of the magnons (•) are de- 
termined by the set (10) of 7,'s and the positions of the 
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FIG. 2: Physical vacuum \G) for a chain of N = 8 spins at 
magnetization M z = 0, 1, . . . ,4. The values of the Bcthe 
quantum numbers Ii are given by the positions of the 
magnons (•). The spinons (O) correspond to vacancies 
in the /^-configurations. 



spinons (O) by the vacancies across the full range of the 
ij's allowed by (S) for class K r states. 

Henceforth we treat \G) as the new physical vacuum. 
At h — (top row in Fig. |2j) it coincides with the spinon 
vacuum, a state with N/2 magnons. At h = hs (bot- 
tom row) it coincides with the magnon vacuum, a state 
containing iV spinons. All states within class K r are gen- 
erated from \G) by rearranging the magnons or (equiva- 
lently) the spinons into all allowed configurations. 

For r = N/2 (top row) and r = (bottom row) the 
state shown is the only possible configuration within class 
K r . In the fourth row, the lone magnon can be moved 
across the array of spinons, generating a branch (one- 
parameter set) of 1-magnon excitations for N — > 00. In 
the second row, the two spinons can be moved indepen- 
dently across the array of magnons, generating a con- 
tinuum (two-parameter set) of 2-spinon excitations for 
N — > 00 with boundaries (Q) as shown in Fig. |l|. The 
center row in Fig. || pertains to the field at half the satu- 
ration magnetization (M z = r = N/4), the case we shall 
investigate extensively for various system sizes. Here \G) 
contains twice as many spinons as it contains magnons. 

The integer m with range < m < M z used in (^) is 
a convenient quantum number for the subdivision of the 
classes K r . Every state of K r at fixed m can then be 
regarded as a scattering state of m pairs of spinon-likc 
quasi-particles. To distinguish them from the spinons, 
we name the new quasi-particles psinons. 

The ground state \G), the only state with m = 0, is the 
psinon vacuum. Here the magnons form a single array 
flanked by two arrays of spinons (see Fig. ||). Relaxing 
the constraint in (^|) from m — to m = 1 yields a two- 
parameter set of states - the 2-psinon excitations. Here 
the array of magnons breaks into three clusters separated 
by the two innermost spinons, which now assume the role 
of psinons. The remaining 2M Z — 2 spinons stay sidelined. 
In the 4-psinon states (m = 2), two additional spinons 
have been mobilized into psinons. By this prescription, 
we can systematically generate sets of 2m-psinon excita- 
tions for < m < M z . 

To illustrate the quasi-particle role of the psinons in 
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FIG. 3: Spectrum of 2-psinon excitations (circles) and 4- 
psinon excitations (squares) for M z — N /4 and N = 16. 
The states marked by full symbols are dynamically dom- 
inant in S zz (q,cu). The spectral ranges of the 2-psinon 
states (thick lines) and 4-psinon states (thin lines) for 
N — > oo are inferred from data for N — 2048. The inset 
shows the psinon energy momentum relation e^(p). 

the class- if,, collective states we have plotted in Fig. H 
energy versus wave number of all 2-psinon states (circles) 
and 4-psinon states (squares) at M z — N/4 for N = 16. 
Also shown are the spectral boundaries of 2-psinon and 
4-psinon states for N — > oo as inferred from solutions of 
(f?J) for N — 2048. The 2-psinon continuum, outlined by 
thick lines, is confined to the interval \q\ < q s , where 

q s EE tt(1 - 2M Z /N) (11) 

denotes the wave number of an incommensurate soft 
mode. The lower 4-psinon spectral boundary is the same 
as the 2-psinon lower boundary but extended periodi- 
cally over the entire Brillouin zone. The upper 4-psinon 
boundary is related to the upper 2-psinon boundary by 
a scale transformation {q — > 2q,tu — > 2u>). 

The relationship between the ranges in (q, w)-space of 
the 2-psinon states and the 4-psinon states does indeed 
reflect the fact that they are scattering states of two 
or four quasi-particles, respectively, of the same species. 
Like the spinon, the psinon is not observable in isolation 
via neutron scattering, but its energy momentum rela- 
tion €ip(p), — 7r/4 < p < 7r/4, can be inferred from the 
data of Fig. || (see inset). 

If there were no psinon interaction, the wave number 
and energy of a 2m-psinon state would be q = Y^n=iPii 
u) = Yli=i e 4>(Pi)- The N — 16 data make it quite clear 
that the finite-size energy correction caused by the psi- 
non interaction is stronger in the 4-psinon states than in 
the 2-psinon states. In both sets of collective states, the 
interaction energy goes to zero as the scattering events 
become less and less frequent in a chain of increasing 



length. However, it takes longer chains for finite-TV 4- 
psinon data to reach comparable convergence toward the 
spectral boundaries predicted for N — > oo, because for 
fixed N, the scattering events between psinons are more 
numerous in a typical 4-psinon state than in a typical 
2-psinon state. 

If instead of the psinon vacuum we had used the spinon 
vacuum as the reference state at M z = N/4, then both 
the 2-psinon states and the 4-psinon states would have 
to be described as scattering states of N/2 spinons. Al- 
though we know the energy-momentum relation of a 
spinon, Eq. (||), it is of little use to determine the spec- 
tral threshold in Fig. ||. Since the 2-psinon and 4-psinon 
states maintain a finite density of spinons in the limit 
N — > oo, the spinon interaction energy remains signifi- 
cant. This problem does not arise at M z — 0. In the 
2-spinon scattering states depicted in Fig. |l|, the spinon 
interaction energy vanishes for N — > oo just as the psi- 
non interaction energy does in the 2-esinon and 4-psinon 
scattering states depicted in Fig. 0£3 

V. DYNAMICALLY RELEVANT EXCITATIONS 

At M z — the spectral weight in the dynamic spin 
structure factor S zz (q,u>) is dominated by the 2-spinon 
excitations.cl Our task here is to determine how the spec- 
tral weight of S zz (q,uj) at M z ^ is distributed among 
the 2m-psinon excitations. In investigating this question, 
we follow the strategy of an older studyO but with vastly 
improved conceptual and numerical tools. 

We begin by exploring, in a chain of N — 16 spins at 
M z /N = |, the transition rates between the ground state 
\G) and all 2m-psinon excitations for m = 0,1,2,3,4. 
The Bethe quantum numbers of the states with m = 0, 1 
are shown in Fig. ^|. The first row represents the psinon 
vacuum with its four magnons sandwiched by two sets of 
four spinons. The two innermost spinons (marked grey) 
become psinons when at least one of them is moved to 
another position. In the rows underneath, the psinons are 
moved systematically across the array of magnons while 
the remaining spinons stay frozen in place. These eight 
configurations describe all 2-psinon states with q > 0. 

The wave numbers, energies, and transition rates of the 
states shown in Fig. Q are listed in Table H. Remarkably, 
almost the entire 2-psinon spectral weight is concentrated 
in the lowest excitation for any given q. The dynamically 
dominant 2-psinon states are marked by solid circles in 
Fig. [| In a macroscopic system, they form the lower 
boundary of the 2-psinon continuum. 

Next we calculate the transition rates | {GIS"^ | A) | 2 for 
the complete set of 4-psinon states. Interestingly, we ob- 
serve that most of the 4-psinon spectral weight is again 
carried by a single branch of excitations. The dynami- 
cally dominant 4-psinon states for N = 16 are shown as 
full squares in Fig. [| For large N they form a branch 
adjacent to the 2-psinon spectral threshold. 
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FIG. 4: Psinon vacuum \G) for N = 16, M z = 4 and 
2-psinon states with q > 0. The 7, values are marked by 
the positions of the magnons (small circles). The spinons 
(large circles) mark Ji-vacancies. A subset of the spinons 
are called psinons (grey circles). 



TABLE I: Ground state \G) and 2-psinon excitations for 
N = 16, M z = 4, and wave numbers q = k — ko > 
(in units of 2ir/N). The ground state has ka = and 
E G = -11.5121346862. 
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An investigation of the remaining 2m-psinon states 
shows that there exists one dynamically dominant branch 
of 2m-psinon excitations for < m < M z . The configu- 
rations of Bethe quantum numbers pertaining to the four 
branches for N = 16, each consisting of N/2 — M z = 4 
states (at q > 0), are shown in Fig. ||. The energies, 
wave numbers, and transition rates of these excitations 
are listed in Table [n]. All other 2m-psinon excitations 
have transition rates that are smaller by at least two or- 
ders of magnitude at q < n/2, and still by more than one 
order of magnitude at q > n/2. 

Inspection of Fig. reveals an interesting pattern, in- 
dicative of the composition of the dynamically relevant 
collective excitations. They form a two-parameter set. 
The two parameters are highlighted by grey circles. Hith- 
erto we have interpreted each group of four configurations 
as a branch of 2m-psinon excitations, which are seem- 
ingly arbitrary one-parameter subsets taken from 2m- 
parameter sets of states. In a macroscopic system, all but 
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: Psinon vacuum \G) for TV 




16, M z = 


= 4 and set 



of iptp* states with < q < it. The U are given by the 
positions of the magnons (small circles) in each row. The 
spinons (large circles) correspond to Ji-vacancies. The 
psinon and the antipsinon (ip*) are marked by a large 
and a small grey circle, respectively. 

the lowest such branches contain a macroscopic number 
of psinons. Hence the range of the dynamically relevant 
excitations in (q, w)-space cannot be inferred from the 
psinon energy-momentum relation alone as was possible 
for the 2-psinon and 4-psinon continua, because the psi- 
non interaction energy will remain non-negligible in most 
of these states for 7Y — > oo, just as the spinon interaction 
energy was non-negligible in the 2-psinon and 4-psinon 
scattering states at M z ^ 0. 

A more natural interpretation of the pattern on dis- 
play in Fig. [5] identifies one of the two parameters as a 
psinon (large grey circle) as before and the other parame- 
ter as a new quasi-particle (small grey circle). The latter 
is represented by a hole in what was one of two spinon 
arrays in the psinon vacuum. Instead of focusing on the 
cascade of psinons (mobile spinons) which this hole has 
knocked out of the vacuum, we focus on the hole itself, 
which has properties commonly attributed to antiparti- 
cles. The psinon (ip) and the antipsinon (ip*) exist in dis- 
junct parts of the psinon vacuum, namely in the magnon 
and spinon arrays, respectively. When they meet at the 
border of the two arrays, they undergo a mutual anni- 
hilation, represented by the step from the second row to 
the top row in Fig. ^. 

We could have interpreted the small grey circle as a 
magnon (spin-1 quasi-particle), but when we do that 
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TABLE II: Ground state and dynamically dominant ex- 
citations for (N = 16, r = 4) among 2m-psinon states 
(to — 0, 1, . . . ,4). The latter form the ipip* continuum in 
the limit N — > oo. The wave numbers q = k — k(j > 
are in units of 2ir/N. 
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FIG. 6: (a) ipip* excitations at M z /N = 1/4 for N = 16 
(circles, squares, diamonds, triangles for m = 1,2,3,4, 
respectively) and N = 256 (dots), (b) Integrated inten- 
sity S zz (q) (inset) and relative ipip* contribution (main 
plot) for N = 12, 16, 20, 24, 28, 32. The lines connect the 
N = 32 data points. 



we must take into account that it then coexists in the 
magnon vacuum with a macroscopic number of fellow 
magnons (small black circles) . From this perspective, the 
collective excitation must be viewed as containing a finite 
density of magnons (for N — > oo), in which the magnon 
interaction remains energetically significant for scatter- 
ing states. The nonzero interaction energy obscures the 
role of individual magnons. 

On the other hand, when the small grey circle is in- 
terpreted as an antipsinon, then it lives in the psinon 
vacuum, i.e. almost in isolation. The only other particle 
present is a psinon (large grey circle). In a macroscopic 
system, the interaction energy in a psinon-antipsinon 
(ipip*) scattering state becomes negligible. Therefore, the 
identity of both quasi-particles is easily recognizable in 
the spectrum. 

The energies versus the wave numbers of the 16 ipip* 
states listed in Table || are shown in Fig. |(|a) as large 
symbols. The four branches from bottom to top pertain 
to to = 1, . . . ,4. Also shown in the same plot are the 
ipip* states for N = 256. The lower boundary of the 
ipip* continuum emerging in the limit N — > oo touches 
down to zero frequency at q = and q = q s = tt/2. 
Between q s and 7T, it rises monotonically and reaches 
the value E — Eq = h. A direct observation of the in- 
commensurate soft mode at q s was made in a neutron 
scattering .experiment on Cu(C6D 5 COO)2-3D 2 (copper 
benzoate) B 

Figure ^(b) shows the relative integrated intensity of 
the ipip* excitations for various N at fixed M z /N = j. At 
1 % Is = t/2, virtually all spectral weight of S zz (q,u) 



originates from ipip* fluctuations. An extrapolation of 
the data points at q — tt/2 suggests that the relative 
ipip* spectral weight is in excess of 93%. 

At g > 5 S the ipip* contribution to the integrated in- 
tensity decreases monotonically but stays dominant over 
more than half the distance to the zone boundary. The 
width of the ipip* continuum vanishes linearly on ap- 
proach of q = 7r, and the relative spectral weight more 
slowly: S zz (q) ~ (it — q) 7 , 7 ~ 0.3. This enhances the ob- 
servability of the ipip* excitations in the narrow energy 
range near the Brillouin zone in spite of the low abso- 
lute intensity. Finite- N data for the integrated intensity 
S zz (q) are shown in the inset to Fig. |^(b). This function 
is peaked at q = q s , where the ipip* spectral weight is 
overwhelmingly predominant. 

When we lower M z , the soft mode at q s moves to the 
right, the number of 2TO-psinon branches that contribute 
to the ipip* continuum shrinks but each branch gains ad- 
ditional states. At M z = 1 we are left with one 2-psinon 
branch extending over the interior of the entire Brillouin 
zone. This branch is equal to the lowest branch of 2- 
spinon states with dispersion ex,(g), Eq. (^). However, 
even for this case the psinon vacuum is different from the 
spinon vacuum. The former is the lowest-energy 2-spinon 
state (with M z — 1), whereas the latter is a state with 
M z = 0. The wave number of the two vacua differ by 
7r. At M z — the ipip* excitations disappear altogether. 
The limit h — ► of the infinite chain is very subtle and 
will be discussed elsewhere.EEl 

When we increase M z toward the saturation value, the 
soft mode moves to the left, and the number of 2TO-psinon 
branches increases, but each branch becomes shorter. At 
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FIG. 7: Energy versus wave number of all ipip* scattering 
states at q > for N = 64 (o) in comparison with the 
corresponding free ipip* states (+). The inset shows the 
energy-momentum relations of the psinon (0 < \p\ < 7r/4) 
and the antipsinon (tt/A < \p\ < 37r/4) as inferred from 
from ipip* data for TV = 2048. 



M z = N/2 — 1, the two-parameter set collapses into a 
one-parameter set consisting of one 2m-psinon state each 
for m = 1,2,... , A/2 — 1. These states are more nat- 
urally interpreted as a branch of 1-magnon excitations 
with dispersion e\(q) = J(l — cosg). Their relative spec- 
tral weight in S zz (q,w) is now 100%, but the absolute 
intensity for q ^ is only of 0(A _1 ). 

To further illustrate the roles of the psinon and the an- 
tipsinon as the relevant quasi-particles in the collective 
excitations dominating the spectral weight in S zz (q,uj), 
we compare in Fig. ^ the energies between the ipip* scat- 
tering states for N = 64 and the corresponding (ficti- 
tious) free ipip* superpositions. The vertical displacement 
of any (o) from the associated (+) reflects the interaction 
energy between the two quasi-particles. This energy ap- 
proaches zero for all states of this class as N — > oo. 

The energy-momentum relations of the two quasi- 
particles can be accurately inferred from N = 2048 data 
for the spectral thresholds of the ipip* states as illustrated 
in the inset to Fig. 0. The psinon dispersion e^,(p) is con- 
fined to the interval at < \p\ < ir/4 (solid line) and the 
antipsinon dispersion e,p* (p) to 7r/4 < \p\ < 37r/4 (dashed 
line). The different ranges of momentum which the two 
quasi-particles are allowed to have correspond to the dif- 
ferent regions in Fig. ^ across which the circles pertaining 
to them can be varied. 

The lower boundary of the ipip* continuum is defined 
by collective states in which one of the two particles has 
zero energy: the psinon for < \q\ < ir/2 and the an- 
tipsinon for 7r/2 < \q\ < tt. The upper boundary consists 
of three distinct segments. 

For < q < 0.3935 the highest-energy ipip* state is 



made up of a zero-energy psinon with momentum p^p = 
— 7r/4 and an antipsinon with momentum p^p- = tt/A + q. 
Here the shape of the continuum boundary is that of the 
psinon dispersion. Likewise, for 37r/4 < q < tt, the states 
along the upper continuum boundary are made up of 
a maximum-energy antipsinon (with momentum p^,- = 
37r/4 and a psinon with momentum p^p = — 37r/4 + q. 
Here the shape of the continuum boundary is that of the 
psinon dispersion. 

When these two delimiting curves are extended into 
the middle segment, 0.3935 < q < 37r/4, they join in a 
cusp singularity at q — ir/2. Here the highest ipip* state 
does not involve any zero-energy quasi-particles. The 
maximum of e^(p^) + e^* (p.,p* ) subject to the constraint 
Pi> + Pip' = q does not occur at the endpoint of any 
quasi-particle dispersion curve. Consequently, the ipip* 
continuum is partially folded about the upper continuum 
boundary along the middle segment. 



VI. 



LINESHAPES 



To calculate the lineshapes relevant for fixed-g scans in 
an inelastic neutron scattering experiment from the spec- 
trum and matrix elements obtained via Bethe ansatz, 
we exploit key properties of transition rates and den- 
sities of states of sets of excitations that form two- 
parameter continua in (g, ui) space for N — > oo. The 
ipip* transition rates (scaled by N) form a continu- 
ous function Mf_f (q,w) for N — > oo. The ipip* den- 
sity of states (scaled by TV^ 1 ) becomes a continuous 
function D^*(q,w) for N — * oo. The ipip* spectral- 
weight distribution is rthen the product S^f (q,u) — 
D^* (q,u)Mff* (q,u)M In the following, we consider 
three wave numbers at M z = N/4. 

At q = 7r/2, the ipip* continuum is gapless and the 
relative ipip* spectral weight in S zz (q, u>) has a maximum. 
The scaled density of ipip* states is generated from N = 
2048 data of the set of points 



2ir/N 



U) U ' + 1 — io v 



(12) 



where v* = m marks the antipsinon quantum number in 
the ipip* continuum and picks the dynamically relevant 
branch from the set of 2m-psinon states. The psinon 
quantum number v is adjusted to keep the wave number 
q of the ipip* state fixed. This choice of labels produces 
an ordered sequence of levels. Starting at id — 0, the 
graph of (tt/2, rises from a nonzero value very 
slowly up to near the upper band edge, where it bends 
into a square-root divergence as shown in Fig. ||(a). The 
divergence is produced by a maximum, of the sequence 
w„. at the fold of the f/;-0*-continuumc3 

In Fig. |8|(b) we show finite- A data at q — tt/2 for the 
scaled transition rates 



Mff {q,^*) ^ N\(G\S z q \ V *)f 



(13) 
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FIG. 8: (a) Density of tpip* states at q — tt/2 evaluated 
via (|l2|) from Bethe ansatz data for N = 2048. (b) Tran- 
sition rates ( |l3| ) between the psinon vacuum and the ipip* 
states at q = tt/2 for N = 12, 16, 20, 24, 28, 32. (c) Line- 
shape at q = 7r/2 of the ipip* contribution to S zz (q,u>). 
All results are for M z = N/4. 



FIG. 9: (a) Density of ipip* states at q = 7r/4,37r/4 eval- 
uated via ( |l2| ) from Bethe ansatz data for N — 2048. (b) 
Transition rates ( [l3|) between the psinon vacuum and 
the iptp* states at q = 7r/4,37r/4 for N = 16,24,32. (c) 
Lineshape at q — tt/4, 37r/4 of the ipip* contribution to 
S zz (q,w). All results are for M z = N/4. 



These data compellingly suggest the existence of a 
smooth function M^f (tt/2,uj) for the tpi/i* transition 
rates in the limit N — > oo, which further highlights the 
physical significance of the psinon and the antipsinon as 
relevant quasi-particles in this situation. The function 
Mf//r (tt/2, ui) is monotonically decreasing with a diver- 
gence at u> — and a cusp singularity at the upper band 
edge wrj ~ 1.679J. 

The product of the transition rate function and the (in- 
terpolated) density of states is shown in Fig. ||(c)B2l The 
curve fitted through the data points represents the ipijj* 
lineshape at q — tt/2 in S zz (q,u>). Its most distinctive 
feature is the double peak due to apparent divergences 
at both band edges. 

The divergence at ui = 0, which is caused by the matrix 
elements, is a power law, oj~ q , with an exponent that is 
exactly kno wn fr pm field theoretic studies of the Heisen- 
berg model.t30'0 For the situation at hand, the value 
is a = 0.4688 .... The divergence at wu is caused by the 
diverging density of states but is weakened if the cusp 
singularity of Mff (tt/2, ui) starts from zero at ui = uijj- 
The expectation is a power-law singularity, (u>u — u>)~P 
with an exponent < (3 < i. 

It is interesting to compare the ipip* transition rate 
function Mff(-rr/2,uj) at M z = N/4 inferred from the 
Bethe ansatz with the 2-spinon transition rate function 
M^ z (tt, ui) at M z = calculated via algebraic analysis.! 
The shape of both functions is similar, but there are some 

(2) 

differences: M zz (tt, ui) has a stronger power-law diver- 
gence at u> = and it approaches zero more rapidly at 
the upper band edge. As a result it produces a monoton- 

(2) 

ically decreasing spectral- weight distribution S zz (tt, u>) 
(see Fig. |l|) notwithstanding the fact that the 2-spinon 



density of states is also a monotonically increasing func- 
tion terminating in a square-root divergence. 

At q — 7r/4 the integrated intensity S zz (q) is only a 
third of what it was at q = tt/2, but spread over a nar- 
rower range of frequencies (see Fig. ^). The bandwidth 
has shrunk to less than a third of the value it had at 
q = tt/2. The relative tpip* contribution to the intensity 
is even larger than at q = tt/2, almost 100%. In this ap- 
plication, the method of analysis is stretched more closely 
to its limits because q — tt/4 exists in fewer manageable 
system sizes. However, the data still make reliable line- 
shape predictions possible. 

The density of states D^f* (tt/4, w), plotted in 
Fig. ^|(a), rises discontinuously from zero to a finite value 
at the spectral threshold, AE ~ 0.379J. From there it 
increases gradually with gradually increasing slopa and 
ends in a cusp singularity at the upper band edge.LJ The 
finite-TV data for the scaled transition rates shown in 
Fig. ||(b) again suggest a smooth w-dependence in the 
form of a monotonically decreasing curve with enhanced 
steepness near both band edges. However, the coun- 
tertrend of the density of states at the upper band edge 
is of sufficient strength to produce a second maximum in 
the line shape again. 

Also shown in Fig. [)] are the corresponding data for the 
ipip* density of states, transition rates, and lineshape at 
q = 3tt/4. Here the relative spectral weight carried by the 
i/jip* excitations is only 83% of the value at q — tt/2, but 
that fraction is concentrated over a frequency band that 
has shrunk to 65% of the width at q = tt/2, while the ab- 
solute intensity remains fairly high (87% of the value at 
q = tt/2). Both quantities, which determine the ipip* line- 
shape, exhibit similar frequency dependences as we have 
already observed for the other two fixed-g scans. The 
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density of states is divergent again at the upper bound- 
ary. The energy gap is now much larger, AE ~ 0.899 J. 
The fact that the lower continuum boundary at q — 37r/4 
coincides with the upper continuum boundary at q = ir/A 
is a consequence of the quasi-particle dispersions as dis- 
cussed previously. 



VII. CONCLUSION 

The spectrum of the completely integrable ID s = \ 
Heisenberg antiferromagnet (|l|) can be generated in more 
than one way from multiple excitations of quasi-particlcs. 
The external magnetic field controls the nature of the 
ground state. In strong fields, it becomes the vacuum of 
magnons and in zero field the vacuum of spinons. The dy- 
namically relevant collective excitations of specific quan- 
tum fluctuations in the two cases are then naturally de- 
scribed as composites of quasi-particles from the respec- 
tive species and are likely to involve only a small number 
of quasi-particles. 

In intermediate magnetic fields, neither the magnons 
nor the spinons provide a useful interpretation of dynam- 
ically relevant collective excitations for the same fluctua- 
tion operators. The ground state itself contains a macro- 
scopic number of quasi-particles from one or the other 
of the two species. However, when it is reconfigured as 
the physical vacuum for psinons and antipsinons, then it 
turns out that the spin fluctuation operator S* induces 
predominantly transitions to ^ip* states, which contain 
just one particle from each kind. 

Like the magnon and the spinon, the psinon and the 
antipsinon are interacting quasi-particles in the Heisen- 
berg model ([!]). In the tpip* scattering states, the inter- 
action energy of the psinon and the antipsinon is of order 
0(N^ 1 ) whereas the interaction energy among magnons 
or spinons is of order 0(1). Hence, for N — > oo, the i/ji/j* 
states join up in (q, oj)-space to form a two-parameter 
continuum whose spectral boundaries and density of 
states are fully determined by the energy-momentum re- 
lations of the psinon and the antipsinon. Moreover, the 
scaled ipip* transition rates converge for N — > oo toward 
a smooth function of q and uj. 

We have exploited these asymptotic quasi-particlc 
properties to extract lineshape information for the dy- 
namic structure factor S zz {q,u), which probes the spin 
fluctuations parallel to the applied magnetic field. The 
same quasi-particles will also play a dominant role in the 
spin fluctuations perpendicular to the field, but here dif- 
ferent combinations of them make up the composition of 
the dynamically relevant collective excitations. In the dy- 
namic spin structure factor S \-(q,tu), for example, the 

spectral weight is almost completely carried by 2-psinon 
excitations.E£ 

In all likelihood, the psinon quasi-particles will also be 
useful for the analysis of thermal spin fluctuations in this 
model system. The peculiar spectral weight distributions 



found in recent complete diagonalization studiesB0 of 
S zz (q,Lu) at h — and T > 0, for example, indicate the 
presence of stringent selection rules between collective 
states coupled by the spin fluctuation operator S*. In 
zero field, psinon vacua are densely spread across the en- 
tire energy range of the model. Each psinon vacuum can 
be used as the reference state of a 2m-psinon expansion 
(^). If there are general selection rules related to psinon 
quasi-particles among transition rates |(A'|Sg|A| 2 within 
a given class K r of Bethe ansatz solutions, they will have 
a strong impact on the spectral weight distribution in 
S zz (q,uj) at all temperatures. 
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